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ABSTRACT 

Aims. Most studies of high-mass star formation focus on massive and/or luminous clumps, but the physical properties of their larger 
scale environment are poorly known. In this work, we aim at characterising the effects of clustered star formation and feedback of 
massive stars on the surrounding medium by studying the distribution of warm gas through mid-/ 12 CO and 13 CO observations. 
Methods. We present APEX 12 CO(6-5), (7-6), 13 CO(6-5), (8-7) and HIFI 13 CO(10-9) maps of the star forming region G327.36-0.6 
with a linear size of ~ 3 pc x 4 pc. We infer the physical properties of the emitting gas on large scales through a local thermodynamic 
equilibrium analysis, while we apply a more sophisticated large velocity gradient approach on selected positions. 
Results. Maps of all lines are dominated in intensity by the photon dominated region around the Hn region G327. 3-0.5. Mid-7 12 CO 
emission is detected over the whole extent of the maps with excitation temperatures ranging from 20 K up to 80 K in the gas around 
the Hn region, and H2 column densities from few 10 21 cnT 2 in the inter-clump gas to 3 x 10 22 cirT 2 towards the hot core G327.3- 
0.6. The warm gas (traced by 12 and 13 CO(6-5) emission) is only a small percentage (~10%) of the total gas in the infrared dark 
cloud, while it reaches values up to ~35% of the total gas in the ring surrounding the Hn region. The 12 CO ladders are qualitatively 
compatible with photon dominated region models for high density gas, but the much weaker than predicted 13 CO emission suggests 
that it comes from a large number of clumps along the line of sight. All lines are detected in the inter-clump gas when averaged over 
a large region with an equivalent radius of 50"(~0.8 pc), implying that the mid-7 12 CO and 13 CO inter-clump emission is due to 
high density components with low filling factor. Finally, the detection of the 13 CO(10-9) line allows to disentangle the effects of gas 
temperature and gas density on the CO emission, which are degenerate in the APEX observations alone. 

Key words. Stars: formation - ISM: Hn regions - ISM: individual objects: G327. 36-0.6 



1. Introduction 

The influence of high-mass stars on the interstellar medium is 
tremendous. During their process of formation, they are s ources 
of powerful, bipolar outflows (e.g., iBeuther et all [2002), their 
strong ultraviolet and far-ultraviolet radiation fields give rise to 
bright Hn and photon dominated regions (PDRs) and during their 
whole lifetime powerful stellar winds interact with the surround- 
ings. Finally, their short life ends in a violent supernova explo- 
sion, injecting heavy elements into the interstellar medium and 
possibly triggering further star formation with the accompany- 
ing shocks. These are also the type of regions that dominate far- 
infrared observations of starburst galaxies. 

Most studies of massive star formation focus on emission 
peaks at infrared or submillimetre wavelengths, which corre- 
spond to peaks in the temperature and/or mass distribution. The 
aim of our work is to characterise the effects of clustered star 
formation and feedback of massive stars on the surrounding 
medium. We have made APEX maps of three cluster-forming 
regions (G327.3-0.6, NGC6334 and W51) in mid-/ 12 CO ((6- 
5) and (7-6)) and 13 CO transitions ((6-5) and (8-7)) in order 
to have a direct measure of the excitation of the warm extended 
inter-clump gas between dense cores in the cluster (see for ex- 



ample iBlitz & Starkl Il986t IStutzki & Giistenl Q990). Our sam- 
ple of sources was chosen among six nearby cluster-forming 
clouds mapped in water and in the 13 CO(10-9) transition as part 
of the Water in Sta r -For ming Regions with Herschel (WISH) 
(Ivan Dishoeck et all I201 lh guaranteed ti me key program (G T- 
KP) for the Herschel Space Observatory (iPilbratt et alll2010l) . 

In this paper, we present the 12 CO and 13 CO maps of 
the star-formin g reg ion G327.3-0.6 at a distance of 3.3 kpc 
dUrquhart et all 1201 ll based on H I absorption). G327. 3-0.6 
is well suited to study cluster-forming clouds because of its 
relatively close distance and because several sources in differ- 
ent evolutionary phase s coexist in a small region, as found by 
IWvrowski et ail (|2006). Our maps (with a line ar extension of 
~ 3 p c x 4 pc) cover the Hn region G327.3-0.5 |Goss & Shaveil 
Il970h associate d with a lumino us PDR, and an infrared dark 
cloud (IRDC) (IWvrowski et all 120061) hosting the bright hot 
core G327.3-0.6 dGibb et all 120001) a nd the extended green ob - 
ject (EGO) candidate G327.30-0.58 (ICvganowski et all 12008). 
EGOs are identified through their extended 4.5 fim emission in 
the Spitzer IRAC2 band, which is beli eved to trace outflows from 
massive young stellar objects (YSOs) (Cvg anowski et all l2008). 
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Fig. 1. Spitzer infrared colour image of the region G327.3-0.6 
with red representing 8.0//m, green 4.5pm and blue 3.6pm. The 
blue contours represent the integrated intensity of the 12 CO(6-5) 
line (thin contours are 15%, 45% and 70% of the peak emission, 
thick contours 30%, 60% and 90% of the peak emission). The 
sources discussed in this paper (the IRDC, the EGO, the hot core 
G327.3-0.6, the SMM6 position and the Hn region G327.3-0.5) 
are also marked with crosses. The white and yellow boxes mark 
the regions mapped with APEX and Herschel, respectively. 



with on-the-fly maps of 200"x240" spaced by 4" in declination 
and right ascension. 

We used the Fast Fourier Transform Spectrometer (FFTS, 
iKlein et aUl2006l) as backend with two units of fixed bandwidth 
of 1.5 GHz and 8192 channels per pixel. We used the two IF 
groups of the FFTS with an offset of +460 MHz between them. 
The original resolution of the dataset is 0.3 km s _1 ; the spectra 
were smoothed to 1 km s for a better signal-to-noise ratio. 

The observations were performed under good weather con- 
ditions with a precipitable water vapour level between 0.5 and 
0.7 mm. Typical single side band system temperatures during 
the observations were around 1600 K and 5200 K, for the low 
and high frequency channel respectively. The conversion from 
antenna temperature units to brightness temperatures was done 
assuming a forward efficiency of 0.95 for all channels, and a 
main beam efficiency of 0.48 for the 12 CO and 13 CO (6-5) ob- 
servations, 0.45 for the 12 CO (7-6) data, and 0.44 for 13 CO(8-7), 
as measured on Jupiter in September 20083 The pointing was 
checked on the continuum of the hot core G327.3-0.6 (aj2ooo = 
15 A 53 m 07*.8,5j2ooo = -54°36'06'/4). The maps were produced 
with the XY_MAP task of CLASS9(fl which convolves the data 
with a Gaussian of one third of the beam: the final angular res- 
olution is 974 for the low frequency data, 8'.' 1 for the high fre- 
quency. 



Table 1. Observational parameters. 



Line 


Frequency 


Telescope 


Beam 


Av 


r.m.s. 




(GHz) 




(") 


(km s -1 ) 


(K) 


CO(6-5) 


691.473 


APEX 


970 


1.0 


1.4 


CO(7-6) 


806.652 


APEX 


777 


1.0 


4.4 


13 CO(6-5) 


661.067 


APEX 


974 


1.0 


1.1 


13 CO(8-7) 


881.273 


APEX 


7/1 


1.0 


2.8 


3 CO(10-9) 


1101.350 


Herschel 


1970 


1.1 


0.1 



This paper is organised as follows: in Sect. |2]we present the 
APEX and Hersche\}\ observations of G327.3-0.6, in Sect. [3] we 
discuss the morphology and kinematics of the 12 CO and 13 CO 
emission, in Sect. [4] we investigate the physical conditions of the 
emitting gas. Finally, in Sect. [5]we discuss our results and com- 
pare to similar observations performed towards low- and high- 
mass star forming regions. Our results are summarised in Sect. [6] 

2. Observations 

2.1. APEX telescope 

The CHAMP+ dKasemann et all [2006I: [Glisten et aUl2008h dual 
colour heterodyne array receiver of 7 pixels per frequency chan- 
nel on the APEX telescopeQ was used in September 2008 to 
simultaneously map the star-forming region G327.3-0.6 in the 
12 CO (6-5) and (7-6) lines and, in a second coverage, the 13 CO 
(6-5) and (8-7) transitions. The region from the hot core in 
G327.3-006 to the Hn region G327.3-0.5 (Fig. [T) was covered 



Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 

2 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
the Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory 



2.2. Herschel Space Observatory 

The 13 CO (10-9) line (see Table [B was mapped (size=210 " x 
270") with the HIFI instrument (ide Graauw et al.L 1201 Oh to- 
wards G327.3-0.6 on February, 18th, 201 1 (observing day (OD) 
645, observing identification number (OBSID) 1342214421. 
The centre of the map is arjaooo = 15 /! 53'"05 s .48, <5j2000 = 
-54°36'06'/2. The observa tions are part of the WISH GT-KP 
(Ivan Dishoeck et all 1201 lh . Data were taken simultaneously in 
H and V polarisations using both the acousto-optical Wide- 
Band Spectrometer (WBS) with 1.1 MHz resolution and the 
correlator-based High-Resolution Spectrometer (HRS) with 250 
kHz nominal resolution. In this paper we present only the WBS 
data. We used the on-the-fly mapping mode with Nyquist sam- 
pling. HIFI receivers are double sideband with a sideband ratio 
close to unity. The double side band system temperatures and to- 
tal integration times are respectively 384 K and 3482 s. The rms 
noise level at 1 km s _1 spectral resolution is ~0.1 K. Calibration 
of th e raw data onto Ta scal e was performed by the in-orbit sys- 
tem ([Roelfsema et aTT. l2.Q12l) : conversion to T m b was done with a 
beam efficiency of 0.74 and a forward efficiency of 0.96. The flux 
scale accuracy is estimated to be around 15% for band 3. Data 
calibration was perf ormed in th e Herschel Interactive Processing 
Environment (HIPE, Ott, 2010) version 6.0. Further analysis was 
done within the CLASS90 package. After inspection, data from 
the two polarisations were averaged together. 

The original angular resolution of the data is 19'.'0. The final 
maps were produced with the XYJVIAP task of CLASS90 and 
have an angular resolution of 2 1'.' 1 . 



3 http://www.mpifr-bonn.mpg.de/div/submmtech/heterodyne/ 
champplus/champmain.html 

4 http://www.iram.fr/IRAMFR/GILDAS/ 
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Fig. 2. Maps of the integrated intensity of the 12 CO(3-2), (6-5) and (7-6) lines in the velocity range Vlsr = [-54, -40] km s _1 
(colour scale). Solid contours are the integrated intensity of the 12 CO(6-5) line from 20% of the peak emission in steps of 10%. In 
each panel, the positions analysed in Sect. 14.21 (the hot core, the IRDC position, and the centre of the Hn region) are marked with 
black triangles. Th e red star labels the position of the EGO. SMM6 (left panel) is one of the submillimetre sources detected by 
iMinieretalJ (120091) . 

i' w on {h. Km s j 

5D 13C 



H a 3&W 




a (JMOO] 



Fig. 3. Maps of the integrated intensity of the 13 CO(6-5), (8-7) and (10-9) transitions in the velocity range Vlsr = [-54, -40] km s _1 
(colour scale). The solid contours in the left p anel represent the LABOCA continuum emission at 870//m from 5% of the peak 
emission in steps of 10% dSchuller et all 120091) . In the middle and right panels, the solid contours are the integrated intensity of the 
13 CO(6-5) from 20% of the peak emission in steps of 10%. In each panel, the positions analysed in Sect. l4.2l (the hot core, the IRDC 
position, and the centre of the Hn region) are marked with black triangles (except in the left panel, where the hot core is shown by a 
white triangle). The red star labels the position of the EGO. 



3. Observational results 

3.1. Morphology 

Figure Q] shows the I2 CO(6-5) integrated intensity map overlaid 
on the composite image of the IRAC Spitzer 3.6, 4.5 and 8.0/im 
bands of the region. The I2 CO emission traces the distribution 
of the 8.0/im emission, but it is also associated with the infrared 
dark cloud found on the east of the hot core. The EGO candi- 



date G327. 30-0.58 identified bv ICveanowski et all d2008l) and 
clearly visible in Fig. Q] is also detected in the I2 CO data as 
secondary peak of emission (Fig. |2j. The map of the integrated 
intensities of the I2 CO(7-6) line is also presented in Fig. [2] to- 
gether with the integrated intensity of the 12 CO(3-2) line from 
IWvrowski et al.l d2006f) . Figure [3] shows the distribution of the 
13 CO(6-5), (8-7) and (10-9) emissions. The accuracy of the rel- 
ative pointings was checked on the hot core G327. 3-0.6. For this 
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Table 2. Coordinates of the main sources in the G327.3-0.6 mas- 
sive star-forming region 



source 


tt J2000 


<5j2000 


SMM6 a 


15:53:00.9 


-54:37:40.0 


hot core b 


15:53:07.8 


-54:37:06.4 


EGO c 


15:53:11.2 


-54:36:48.0 


Hn d 


15:53:03.0 


-54:35:25.6 



a | MinieretalJd2009h 
b Wvrowski et alJ d2008|) 
b Wvrowski et"ai] J2OO8I) 
c Icyganowski et alJ d2008) 

d peak of the centimetre continuum emission from ATCA archival 
data at 2.3 GHz, project number C772 



purpose, we derived integrated intensity maps of lines detected 
only towards this position, which are close in frequency to the 
12 CO(6-5), I3 CO(6-5) and 13 CO(8-7) transitions, and therefore 
were observed simultaneously to the current dataset. From these 
data, we infer a position for the hot c ore in agreement with i nter- 
ferometric measurements at 3 mm (IWvrowski et all 120081 and 
Table © within ~1'.'5. 

All observed I2 CO transitions trace the Hn region G327.3- 
0.5 as well as the infrared dark cloud which hosts the hot core 
G327.3-0.6. Moreover, the 12 CO(6-5), (7-6) and 13 CO(6-5) 
lines show extended emission along a ridge running approxi- 
mately N-S that matches very well with the distribution of the 
CO(3-2) transition. The hourglass shape hole to the west of the 
Hn region G 327. 3-0.5 where t he 12 C O(3-2) emission is strongly 
reduced (see Wvro wski et all 120061) is seen also in the I2 CO(6- 
5) and (7-6) lines which, although much weaker than in the rest 
of the map, are still detected at this position. 

All transitions peak towards the Hn region G327.3-0.5 
where the main isotopologue lines have intensities up to 60- 
65 K. The integrated intensities of the CO isotopologues show a 
distribution along a ri ng-like structure aroun d the peak of the cm 
continuum emission dGoss & ShaverL Il970h . The centre of the 
ring also coincides with the massive young stellar object n umber 
87 identified in the near-infrared by Moise s~t alJ d201 lb . Since 
the ring is detected also in high-/ transitions of 13 CO, it is plau- 
sible that this morphology is true and not due to optical depth 
effects. This structure likely coincides with the limb brightening 
of the hot surface of a PDR around G327.3-0.5 and could trace 
an expanding shell. We will investigate this scenario in Sect. 13.21 

The hot core G327. 3-0.6 shows up as a secondary peak in the 
integrated intensity maps of the 13 CO transitions, while the main 
CO isotopologue peaks to its north-west, probably because of 
optical depth effects. Strong self-absorption profiles are indeed 
detected in all 12 CO lines towards the hot core (see Sect. 13.21 1. 
The submillimetre source SMM6 (see n in the continuum emis- 
sion at 450yum by Minieret al., 2009) is detected as a peak of 
emission in all integrated intensity maps of 12 CO and in I3 CO(6- 
5), although at the edge of the mapped region. The other submil- 
limetre sources are also marked in Figs. [2] The EGO candidate 
G327.30-0.58 is also detected in the I3 CO(6-5) map (Fig. |3). 
The 13 CO(6-5) traces the whole IRDC and not only the active 
site of star formation where the EGO is detected. The contin- 
uum emission due to dust (seen for example at 870/mi in Fig. [3) 
follows the distribution of the 13 CO lines. 

In Fig. 3] we show the ratio of the integrated intensity of 
the 12 CO(6-5) transition (convolved to the 18" resolution of the 
12 CO(3-2) data) to that of the 12 CO(3-2) line. This ratio ranges 




a [J2D00] 



Fig. 4. Distribution of the line ratio of the 12 CO(6-5) transition 
to the 12 CO(3-2) line. Solid contours show the 12 CO(6-5) inte- 
grated intensity in the velocity range Vlsr = [-54, -40] km s _1 
from 20% of the peak value in steps of 10%. The black triangles 
are as in Fig. [2] 



between 0.3 and 1.8; it has values slightly larger than one to- 
wards the Hn region (1.2 at its centre), while it is about unity 
towards the hot core. The peak is found south-west of the Hn 
region G327. 3-0.5, where both lines are detected with a high 
confidence level. However, these results could be biased by the 
strong self absorption in both I2 CO lines (see Sect. 13 .21 . For this 
reason, we computed the ratio between the two transitions in 
four velocity ranges to cross-check the results of Fig. [4] The in- 
ferred values, however, do not change significantly. 

3.2. Line profiles and velocity field 

The widespread 12 CO(6-5) emission shows line profiles with 
a typical width of ~8 km s _1 in the gas between G327. 3-0.5 
and the infrared dark cloud. Broader profiles are detected in the 
infrared dark cloud and in the northern part of the Hn region 
G327.3-0.5. Figure [5] shows the distribution of the line width 
of the 12 CO(6-5) transition: the 12 CO(6-5) line width follows 
an arc-like structure that connects the Hn region G327. 3-0.5 to 
the infrared dark cloud where the hot core is. Interestingly, the 
same morphology is seen in the LABOCA map of the region 
dSchuller et alll2009h . Line widths are similar for all 12 CO lines, 
while they are consistently narrower in the 13 CO transitions. 

Representative spectra of all 12 CO transitions analysed in 
this study are presented in Fig.|6]towards the hot core, the IRDC 
position ((30", 30") from the centre of the APEX maps, see 
Sect. l4.2l and Figs. [2](3j and the peak of the cm continuum emis- 
sion in G327.3-0.5. Spectra of the 13 CO transitions are shown in 
Fig lZ] Red- and blue-shifted wings are detected in the 12 CO lines 
in a velocity range between -71 and -24 km s _I (in I2 CO(6-5)) 
towards G327. 3-0.6 probably due to outflow motions. However, 
no sign of bipolar outflows is found when inspecting the in- 
tegrated intensity maps of the blue- and red-shifted wings nor 
in position- velocity diagrams (Fig. [HJ. Moreover, very similar 
broad lines are detected along the whole extent of the infrared 
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Fig. 5. Distribution of the second moment of the 12 CO(6-5) line. 
The black solid contours represent the LABOCA continuum 
emission at 870/mi from 5% of the peak emission in steps of 
5%. The peak of the continuum emission corresponds to the hot 
core position. 



dark cloud, as shown in the top panel of Fig. [8] At the IRDC po- 
sition, the wings in 12 CO(6-5) range from -60 to -31 km s . All 
main isotopologue transitions analysed in this paper are affected 
by self-absorption (see Fig. [6] for reference spectra towards the 
hot core, the IRDC position and the Hn region); moreover, even 
the 13 CO(6-5) line shows weak evidence of self-absorbed pro- 
file towards the hot core. Figure [9] shows the 12 CO(7-6) spectra 
overlaid on the continuum emission at 870pm: the self-reversed 
profile is spread over a large area and seems to follow the ther- 
mal dust continuum emission. 

Finally, the velocity field of the 12 CO transitions may help 
us to understand the nature of the ring detected towards the Hn 
region G327.3-0.5. We therefore used the task KSHELL buil t 
in the visualisation software package KARMA (Gooch, 1996). 
KSHELL computes an average brightness temperature on annuli 
about a user defined centre. A spherically symmetric expand- 
ing shell will look like a half ellipse in a f r-v) diagram with the 
axis in the v direction twice the expansion velocity. Figure \W\ 
shows the resulting (r-v) diagram obtained with the 12 CO(6-5) 
data cube using the peak of the cm continuum emission as cen- 
tre. The emission does not follow a perfect spherical shell. This 
is likely due to inhomogeneities in the distribution of the gas, 
as already seen in Fig. QT| where the distribution of the optical 
depth of 13 CO is not symmetric. 



4. Physical conditions of the warm gas 

4.1. LIE analysis 

From the line ratio of the 12 CO(6-5) to 13 CO(6-5) transitions 
we can derive the optical depth of the 12 CO(6-5) line emission, 
which can be then used to infer the excitation temperature of the 
line and the column density of 12 CO in the region. 





20 




-80 -60 -40 -20 
Velocity i.kix/51 

Fig. 6. Mean spectra over a beam of 18" (to match the resolu- 
tion of the 12 CO(3-2) data) of the I2 CO isotopologue transitions 
analysed in the paper. The top panel shows spectra at the centre 
of G327. 3-0.5, the middle panel spectra at the IRDC position, 
the bottom panel spectra at the hot core position. 



The line intensity in a given velocity channel of a given tran- 
sition is 



T L = n x [F v (T ex ) - F v (T cbg )] x (1 - e~ T ") 



(1) 



where rj is the beam filling factor (assumed to be 1 in the follow- 
ing analysis), F v - hv/k x [exp(hv/kT) - r c b g = 2.7 K, and 
t v is the optical depth. Under the local thermodynamic equilib- 
rium (LTE) assumption, T ex is assumed to be equal to the kinetic 
temperature of the gas and equal for all transitions. In the fol- 
lowing analysis, we study the peak intensities of the I2 CO(6-5) 
and 13 CO(6-5) lines, and include only the cosmic background as 
background radiation and neglect, for example, any contribution 
from infrared dust emission since we do not have any map of the 
distribution of the dust temperature. This most likely affects only 
our estimates at the hot core position and possibly towards the 
Hn region G327. 3-0.5 where SABOCA continuum emission at 
350/im is also detected (Wyrowski et al., in prep.). Fo r an appro- 
priate analysis of the emission from the hot core, see lRolffs et al.l 
(1201 lh . 

Assuming that the I2 CO(6-5) emission is optically thick and 
that the 12 CO(6-5) and 13 CO(6-5) lines have the same excita- 
tion temperatures, the optical depth of the I3 CO(6-5) transition, 



Ti3 C0 , IS 



Ti3 C o = -ln\l- 



T L ( l3 CO) 
T L ( n CO) 



(2) 



The optical depth of the 12 CO(6-5) transition can then be 
obtained by multiplying for the abundance of 12 CO relative to 
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Fig. 7. Mean spectra over a beam of 21'.'1 (to match the resolu- 
tion of the 13 CO(10-9) data) of the 13 CO isotopologue transi- 
tions. The selected positions are those discussed in Sect. 14.21 



13 CO, Xi2 CO /i3co ~ 60 dWilson & Roodl 11994 . From the opti- 
cal depth of the 12 CO(6-5) line, one can also derive its excita- 
tion temperature using Eq. Q] Figure QT| shows the distribution 
of the optical depth of the 13 CO(6-5) line and of the excitation 
temperature of 12 CO(6-5). The 13 CO(6-5) emission is moder- 
ately optically thick (0.6-0.7) at the Hii region and at the in- 
frared dark cloud, while it reaches values of ~ 1 .2 at the hot core 
position and in a small part of ring around the Hii region. The 
map distribution of the excitation temperature of the 12 CO(6-5) 
line is shown in the bottom panel of Fig. QT| The map is domi- 
nated by the Hii region, where T ex reaches values of 80 K in the 
ring around the Hii region and then decreases with increasing 
distance from it. The hot core and the rest of the infrared dark 
cloud have values around 30-35 K. The excitation temperature 
increases to the south west of the hot core, in a region where 
there is also 8/im emission, and to the north-east of the Hii along 
a layer of gas also visible in the 12 CO(6-5) integrated intensity 
map (see Fig. O, but more prominent in the T ex map and in the 
8/im emission map (see Fig.Q]). 

From the optical depth and the excitation temperature of 
the 12 CO(6-5) line, we derived the H2 column density assum- 
ing a relative abu ndance of 12 CO relative to H2 of 2.7 x 10 4 
(lLacv et allfl99l . Results are shown in Fig. [12] The largest col- 
umn density is found towards the hot core (~3 x 10 22 crrT 2 in 
the 9'.' 4 beam of the 13 CO(6-5) data) and decreases along the in- 
frared dark cloud with a distribution similar to that of the 870yum 




Fig. 8. Top panel: Colour scale and contours show the P-V di- 
agram of the CO(6-5) transition computed along the cut indi- 
cated by the white arrow in the bottom panel. Offset positions 
increase along the direction of the arrow shown in the bottom 
panel. Bottom panel: distribution of the integrated intensity of 
the 13 CO(6-5) transition towards the hot core G327. 3-0.6. Solid 
contours show the continuum emission at 350/mi (Wyrowski et 
al., in prep.) from 3<x in steps of 10cr (<x ~ 3 Jy/beam). Symbols 
are as in Fig. [2] 



continuum emission. Three peaks around 10 22 cm 2 are found in 
the Hn region. 

We cross-checked our results by computing the H2 column 
density under the assumption that the 13 CO emission is opti- 
cally thin. The results are consistent with those presented in 
Fig. [TJl the largest differences (of the order of 30%) are found 
towards those positions where the optical depth of the 13 CO(6- 
5) transition (Fig. ITTb exceeds ~0.7. The assumption of opti- 
cally thin emission for 13 CO may be particularly useful for the 
inter-clump medium (arbitrarily defined as the region in the map 
where the I3 CO lines are not detected on individual spectra), to- 
wards which we infer H2 column densities of ~2 x 10 21 cm -2 
corresponding to a 13 CO column density of ~10 16 crrT 2 (see also 
Sect. ED- 

We also computed column densities and rotational tempera- 
tures in the region using the rotational diagram technique applied 
to the 13 CO data. We did not include the I2 CO lines in the analy- 
sis because of their complex line profiles and high optical depths. 
Given the optical depth previously derived for the 13 CO(6-5) 
line, we did not apply any correction due to optical depth effects 
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Fig. 9. Map of the 12 CO(7-6) line overlaid on the continuum 
emission at 870jum from LABOCA. The velocity axis ranges 
from -60 to -25 km s _1 , the temperature axis from -1 to 15 K. 
The I2 CO(7-6) data were smoothed to a resolution of 18" to 
match the resolution of the I2 CO(3-2) data and of the LABOCA 
emission. The centre of the map is that of the APEX data (see 
Sect. 12. 11 1. The triangles and the green star are as in Fig. [2] 



to the 13 CO data. Results are consistent with the estimates based 
on equation Q] The main differences between the two analyses 
are found towards the hot core, where the rotational temperature 
is higher than the excitation derived with equation[T](r iot ~70 K 
and T ex ~32 K) and the H2 column density lower (10 22 cm -2 
versus 3 x 10 22 cm -2 obtained with the first method). The dif- 
ferences between the two methods are likely influenced by the 
self-absorption profile detected in the 12 CO(6-5) line which re- 
sults in lower line intensities towards the regions of large column 
densities. In particular, while the first approach overestimates the 
optical depths (and hence column density) because of the self ab- 
sorption, the column density derived with the rotational diagram 
analysis is likely underestimated towards the hot core because of 
the optically thin emission assumption, and should be corrected 
by a factor r LTE /(l - exp(r LTE )-) -1.7. 

Finally, we computed the ratio between the amount of warm 
gas (traced by 12 CO and 13 CO(6-5) and shown in Fig. [T2l 
and the total amount of gas (traced by the continuum emis- 
sion at 870ji/m) as derived assuming a dust temperature equal 
to the excitation temperature of 12 CO(6- 5) and a dust opacity of 
0.0182 cm g 1 feauffman rjetaTl 12008). The results are shown 
in Fig. [13] the warm gas is only a small percentage (~10%) of 
the total gas in the infrared dark cloud, while it reaches values up 
to ~35% of the total gas in the ring surrounding the Hn region. 

4.2. n CO ladder 

Since the G327.3-0.6 region was mapped in three different tran- 
sitions of the 13 CO molecule, we can perform a multi-line anal- 
ysis towards selected positions and infer the parameter of the 
gas. The advantage of 13 CO compared to the main isotopo- 
logue is in the lower opacities of the lines and in the less com- 
plex li ne profiles. For this ana lysis, we used the RADEX pro- 
gram dvan der Tak et al.L 120071) with expanding sphere geome- 
try. The molecular dataset comes from the LAMDA database 




10 20 30 40 
Radius (arcsec) 

Fig. 10. (V-v) dia gram of the Hn region G327.3— 0.5 obtained 
from the 12 CO(6-5) data cube. The radius axis is the distance 
to the shell expansion centre, chosen to be the peak of the cm 
continuum emission. The half ellipse represents an ideal shell in 
( r-v) diagram with an expansion velocity of 5 km s _I . 

Table 3. Line parameters of the 13 CO lines 



Position 


VLSR 

(km s" 1 ) 


FWHM 

(km s" 1 ) 


J T mb 6v 
(Kkms- 1 ) 






1J CO(6-5) 




hot core 


-44.74 ± 0.03 


7.23 ± 0.06 


107.8 ± 0.8 


IRDC 


-45.99 ± 0.02 


5.44 ± 0.05 


67.3 ± 0.8 


Hn 


-48.97 + 0.01 


8.40 ± 0.02 


202.2 ± 0.5 






13 CO(8-7) 




hot core 


-44.63 ± 0.06 


6.5 ± 0.2 


43.5 ±0.9 


IRDC 


-46.5 ± 0.2 


4.4 ± 0.6 


16 ± 2 


Hn 


-48.83 ± 0.05 


8.2 + 0.1 


1 16 ± 1 






13 CO(10-9) 




hot core 


-44.28 ± 0.07 


6.4 + 0.2 


24.4 ± 0.5 


IRDC 


-46.8 ±0.1 


3.7 ±0.3 


4.7 ±0.3 


Hn 


-49.56 ± 0.04 


7.67 ± 0.09 


52.4 ±0.5 



dSchoier et al.Ll2005l) and includes collisional rates adapted from 
lYang et all d2010l) . We ran models with temperatures from 20 to 
200 K, densities in the range 10 4 - 5 x 10 7 cm -3 , and 13 CO 
column densities between 10 14 and 10 19 crrT 2 . 

All data were smoothed to the resolution of the I3 CO(10-9) 
map. We selected three positions for the analysis: the hot core, 
the IRDC position ((30", 30") from the centre of the APEX 
maps), and the centre of the Hn region. The IRDC position was 
selected to be a position associated with high column density in 
the infrared dark cloud (see Figs. [2] [3land[T2l but without IR 
emissio n. However, it is only 10 " to the north of the EGO can- 
didate dCvganowski et all I2008L see Figs. |2]|3]l, and therefore, 
given the beam of the observations, contamination from the em- 
bedded YSO may still be possible. Table |3]reports the measured 
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Fig. 11. Distribution of the optical depth of the 13 CO(6-5) line 
(ti3 CO , top panel) and of the excitation temperature of 12 CO(6-5) 
transition (T ex co, bottom panel). Black contours are the 13 CO(6- 
5) integrated intensity as in Fig. [3] 



line parameters of the 13 CO transitions obtained with Gaussian 
fits; based on these values, we adopt line widths of 6, 3 and 7, 
for the hot core, the IRDC a nd the Hn respectively , in a gree- 
ment with values reported by ISan Jose-Garcia et alJ (j 2012) for 
the 13 CO(10-9) line towards a sample of intermediate- and high- 
mass sources. The spectra are shown in Fig. [7] 

The results of the RADEX analysis are listed in Tableland 
shown in Fig.Qj] For the IRDC position, the 13 CO(10-9) line 
intensity is not well fitted by our one-temperature model. This 
likely reflects the fact that the 13 CO(10-9) spectrum is domi- 
nated by the embedded source while the other two lines sample 
colder gas in the envelope. We considered a 20% calibration er- 
ror for the 13 CO(6-5) a nd (8-7) observations a nd a 15% error 
for the 13 CO(10-9) data. I Wvrowski et a l. (2006) mapped the re- 
gion with the APEX telescope in the C ls O(3-2) line. Therefore, 
since no observations were performed in the 13 CO(3-2) tran- 
sition, we included the C 18 0(3-2) data in Fig. [14] Note how- 
ever, that the C ls O(3-2) fluxes are not included in the fitting 
procedure, but that they are simply used to cross-check results. 
Since the C 18 0(3-2) flux corresponds to a lower limit to the 
flux of 13 CO(3-2) line, we also plotted the C 18 0(3-2) flux cor- 
rected for the abundanc e ratio of 13 CO to C ls O, Xi3 CO /c 18 o ~ 8 
(IWilson & RoodL [r994). This value is likely an upper limit to the 
flux of the 13 CO(3-2) line due to opacity effects. 




50 -50 -100 

Ac ["] 



Fig. 12. Distribution of the H2 column density in the G327. 3-0.6 
star-forming region based on equation[TJ Black contours are the 
13 CO(6-5) integrated intensity as in Fig. [3] 




5d -SSt -100 
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Fig. 13. Distribution of the ratio between the column density of 
warm gas (traced by I2 CO and 13 CO(6-5)) and the total H2 col- 
umn density (traced by the continuum emission at 870/mi) in the 
G327.3-0.6 star-forming region. Black contours are the 13 CO(6- 
5) integrated intensity as in Fig. [3] 



An example of the^- 2 distribution projected on to the T — n 
plane is shown in Fig.[T5]for the HII position, where the reduced 
X 2 at the best fit position is 3. Figure [TBI shows the typical in- 
verse n - T relationship often seen in x 2 distributions and due 
to the fact that density and temperature are, in the case of the 
CO molecule, not indep endent parameters (see Appendix C of 
Ivan der Tak eTalL 120071) . 
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Fig. 14. Distribution of the 13 CO peak line intensities. Full black 
triangles correspond to 13 CO(6-5), (8-7) and (10-9) observed 
intensities. The circle represents the observed C I8 0(3-2) flux, 
the empty black triangle the flux of the C 18 0(3-2) line multi- 
plied by Xi3co/c 18 o ~ 8- The red triangles are the best model fit 
results. The error bars include only calibration uncertainties. The 
dashed lines represent the best fit 13 CO ladder. 



We note here that the detection of the 13 CO(10-9) line breaks 
the degeneracy between density and temperature typical of 1 2 CO 
analyses (e.g., Kra mer et al.L 120041 Ivan der Tak et an,l2007l) and 
help to give stronger constraints: indeed, with the exception of 
the hot core, at least the temperature of the gas is well determined 
at all positions. 

The results of the RADEX analysis confirm that the LTE as- 
sumption used in Sect. 14. H is reasonable since the inferred densi- 
ties are much larger than the critical densities (a few 10 4 cm" 3 for 
all three analysed transitions). Assuming Xi2 CO /n co ~ 60 and an 
abundance of 12 CO relative to H 2 of 2.7 x 10" 4 , the derived 13 CO 
column densities listed in Table|4]correspond to H2 column den- 
sities of some 10 22 cm" 2 for the hot core and the IRDC position, 
and of 2 x 10 22 cm -2 for the Hn position. The derived column 
densities and temperatures are in agreement with those derived 
in Sect. 14. II through rotational diagrams of the 13 CO emission. 
The 13 CO(6-5) optical depths are also in agreement with those 
estimated in Sect. @7T] through the 13 CO and 12 CO(6-5) line ra- 
tio, with the exception of the hot core position: tlte ~1.2 and 
tradex ~0.3 for the hot core, 0.6 and 0.5 for the Hn region, and 
0.6 and 0.7 for the IRDC position. 

We finally stress that the results obtained with the LTE analy- 
sis (Sect. l4~TT ) and with the RADEX code (this section) are based 
on the assumptions that 1) all lines have a beam filling factor of 
one, and 2) that the emitting gas is homogeneous, whereas self- 
absorption profiles in the 12 CO lines indicate an excitation gradi- 
ent along the line of sights. For optically thin lines ( 13 CO in the 
current case), a beam filling factor less than one (but equal for all 
transitions) would mostly affect the column density and result in 




Fig. 15. Projection of the 3-dimensional (T-n-N) distribution of 
the^ 2 on the T — n plane for the Hn position. The contours show 
the 1, 2 and 3<x confidence levels for two degrees of freedom. 
The triangle marks the best fit position. 

Table 4. Best fit parameters of the 13 CO line modelling 



Position 


T 


»H 2 


Mia, 




(K) 


(cm" 3 ) 


(cm" 2 ) 


hot core 


70(> 60) 


1 x 10'(> 10 4 ) 


5 x 10 lb 


IRDC 


45 +20 


1 x 10 6 (> 10 5 ) 


5 x 10 16 


Hn 


75™ 


5 x 10 6 (> 4x 10 4 ) 


1 x 10 17 



The errors represent the 3cr confidence levels in the temperature-density 
plane or 3<x lower (shown in brackets) limit when no stronger con- 
straints can be inferred. 



larger values of N; for optically thick lines, a smaller value of r\ 
would imply larger values of density and/or temperature. 

The uncertainties on the derived parameters due to the as- 
sumption of a homogeneous medium are of less immediate in- 
terpretation, and more complex r adiative transfer codes (e.g., 
iHogerheiide & van der TakL 12000) should be used to reproduce 
the observed line velocity profiles. 



5. Discussion 

5.1. Total 12 CO and 13 CO emission 

Figure[l6]shows the 12 CO and I3 CO ladders obtained by averag- 
ing the emission of the different observed transitions, smoothed 
to the resolution of the HIFI data, over four regions: the total 
map, the Hn region G327.3-0.5, the IRDC hosting the hot core 
(the selected region does include the hot core), and finally the 
inter-clump gas, which was defined as the region in the map 
where the 13 CO lines are not detected on individual spectra. This 
region has an equivalent radius of 50"(corresponding to ~0.8 pc 
at the distance of the source). Examples of 12 CO line profiles to- 
wards the inter-clump gas are shown in Fig. [17] The C 18 0(3-2) 
cannot be used in this analysis because the observations cover a 
much smaller region than that mapped in 13 CO. The total mass 
of the mapped region can be computed using the excitation tem- 
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Fig. 16. CO (left) and 13 CO (right) ladders for the whole mapped region (black), the Hn region G327. 3-0.5 (red), the IRDC 
including the hot core (blue) and the inter-clump gas (magenta). Error bars inclu de only calibration un certainties. In both panels, 
the dashed and dotted curves represent the predicted intensities for model B from lKoester et all ( 1 19941) for a density of 10 7 cm 3 , 
incident UV fields of 10 3 and 10 4 relative to the average interstellar field 

, a visual extinction of 10, and a Doppler broadening of 3 (dashed curve) and 1 km s -1 (dotted curve). 
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Fig. 17. CO(6-5) (solid line) and 12 CO(7-6) (dashed line) spec- 
tra towards some positions in the inter-clump gas. The offset po- 
sition from the centre of the APEX 12 CO maps (Sect. 12. Il l is 
shown for each spectrum in the top right corner. 



perature and H2 column density distributions shown in Figs.fTTI 
and[T2] This corresponds to ~ 700 M . 

The four regions have similar 12 CO and 13 CO ladders. In 
order to correct for self-absorption, integrated fluxes were ob- 
tained for each of the four regions from line fitting of the CO 
spectra with one Gaussian component. While this works fine for 
the spectra of the IRDC hosting the hot core and for the inter- 



clump gas, the line profiles of the total map and of Hn region are 
red-skewed and therefore the fluxes derived with this method 
are likely underestimated. For the main isotopologue, the flux of 
the (7-6) and (6-5) transitions is very similar, although for the 
IRDC+HC and the inter-clump the flux of the (7-6) line is lower 
than that of the (6-5) line. On the other hand, in the 13 CO lad- 
der the peak flux decreases with increasing energy level. For all 
transitions presented in this paper, the spectra are dominated in 
intensity by the Hn region, whose flux is of the order of 60% of 
the total flux for the main isotopologue lines, and ranging from 
~ 49% of the total flux in the (6-5) line to ~ 85% in the (10-9) 
transition for 13 CO (see Table |5j. The main difference between 
the 12 CO spectra and those of 13 CO lies in the emission from 
the inter-clump gas: for all 12 CO lines, the intensity is relatively 
strong, 10% of the flux from the whole region. On the other 
hand, the flux of the 13 CO lines coming from the inter-clump 
gas decreases with increasing J, from ~ 7% of the total flux for 
J=6 to ~ 2% for J=10. The general behaviour of the 12 CO and 
13 CO ladders i s qual itatively compatible with PDR models from 
iKoester et ail d 19941) for high density (10 6 - 10 7 cm 3 ) clouds 
illuminated on one side by a UV radiation field (their model 
B). In Fig. [16] we show the predicted CO and 13 CO line in- 
tensities for models with a density of 10 7 cirT 3 , incident UV 
fields wit h strength 10 3 and 10 4 relative to the average interstel- 
lar field (DraineT. ri978l) . a visual extinction of 10, and a Doppler 
broadening of 3 and 1 km s . High densities (« > 10 6 cm -3 , 
in agreement with our results from Sect. 14.21 1 are needed to lo- 
cate the peak of the CO ladd ers at mid- , high-/ transitions (see 
Figs. 9-10, 12-13 of lKoesteretalJl994h . Stronger UV radiation 
fields also shift the peak of the CO ladders to higher J transi- 
tions than observed. In Fig. [16] the 12 CO model intensities are 
corrected by a factor 0.2 and 0.25 to correct for different line- 
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Table 5. Percentage of integrated fluxes from the Hn region 
G327.3-0.5 and from the inter-clump gas respect to the total in- 
tegrated flux in the map. 



Line 


Hn 


Inter 


CO(3-2) 


56% 


10% 


CO(6-5) 


58% 


10% 


CO(7-6) 


60% 


10% 


13 CO(6-5) 


49% 


7% 


13 CO(8-7) 


67% 


5% 


I3 CO(10-9) 


85% 


2% 



widths between the model and the observations, and possibly 
for beam dilution effects. On the other hand, the 13 CO results of 
Fig.[l6]are not scaled down b y any factor as they are far too weak 
to match the observations. iKoester et al.l (1 19941) already noticed 
that the predicted line intensities of mid-/ 13 CO transitions in 
their models are much weaker than observed in star-forming re- 
gions. They proposed that mid-7 13 CO emission comes from a 
large number of filamentary structures, or clumps, along the line 
of sight. In this way, the modelled line intensity of mid-7 13 CO 
lines would increase significantly as the lines are optically thin, 
while it would not change for optically thick transitions. 

Rotational diagrams of the 13 CO emission applied to the 
spectra of the Hn region, of the IRDC and of the inter-clump 
gas infer rotational temperatures of 66 K, 47 K and 44 K, re- 
spectively, and I3 CO column densities of 6 x 10 15 cm -2 for the 
Hn region and the IRDC, and of 2xl0 15 cirT 2 for the inter-clump 
gas. Since the inter-clump gas has physical parameters very sim- 
ilar to those of the IRDC region, but a much lower column den- 
sity, we suggest that it is composed of high-density clumps with 
low fillin g factors. Th i s is again in agreement w ith PDR mod- 
els (e.g jKoester et aUll994flCubick et alll2008l) which predict 
strong emission at mid-7 13 CO and high-7 12 CO lines in the case 
of small, low ma s s, hi gh density clumps. 

ICubick et al] d2008l) suggested that the COBE 12 CO ladder 
of the Milky Way can be reproduced by a clumpy PDR model, 
and that the bulk of the Galactic FIR line emission comes from 
PDRs around the Galactic population of massive stars. Our ob- 
servations seem to confirm this result, since the CO emission 
of the G327.3-0.6 region is dominated by the PDR around the 
Hn region. Our resu lts are also consisten t with t he findings from 
iDavies et ail (120111) and iMottram et"all (120111) . These authors 
studied the propertie s of massive YSOs and compact Hn regions 
in the RMS survey dHoare et al ., 2005), and found that there is 
no significant population of massive YSOs above ~ 10 5 L Q , 
while compact Hn regions are detected up to ~ 10 6 L . Since 
high-7 CO lines are among the most important cooling lines in 
PDR, they reflect the luminosity of their heating sources: if the 
luminosity distribution of massive stars in the Galaxy is domi- 
nated by Hn regions and not by younger massive stars, then we 
also expect that the CO distribution follows the same rule. 



5.2. Comparisons with other star forming regions 

Large-scale mapping of some low- and high-mass star forming 
regions was performed in several 12 CO transitions. However, 
given the different critical densities of the lines, we prefer to 
compare our results with studies carried on with J> 3 CO lines. 
Excitation temperatures around 8-30 K are found in extended 
diffuse emission and towards the brightest positions, respec- 
tively, in low-mass star forming regions (e.g jDavis et afll2010t 
ICurtis et all 12010). while massive star-forming regions are usu- 



ally warmer as confirmed by our findings (see also Wilson et al 
1 20011: Uakob et all 120071: Ikirk et all l2010L IWilson et all 1201 1 
Pen g et all 120121) . The CO column densities (10 17 - 10 18 cm' 2 
from the LTE analysis) found in G327 are also compara- 
ble w ith values obtained in other massive star-forming re- 
gions (IWilson et all 120011: Uakob eTall 120071: IWhite et al., 2010; 
IWilson et all 1201 lh . 

Comparison with 12 CO large maps of other high-mass star 
forming regions would be important to verify whether our re- 
sult that the 12 CO distribution is dominated by the Hn region 
is a common feature or not. However, most studies do not 
co ver different evol utionary phases as in our case. For OMC- 
1, iPeng et all d2012l) confirmed that the peak of the integrated 
intensity of several CO isotopologue lines is close to the Orion- 
KL hot core (although Orion-south and the Orion Bar PDR are 
also very prominent). One should notice however, that Orion is 
roughly six times closer to the Sun than G327.3-0.6, and there- 
fore the Orion Bar and Orion-KL would be much closer on sky 
(~ 30") if one would place them at the distance of G327.3- 
0.6. Moreover, Orion-KL likely repres ents a special case since 
it hosts a very powerful outflow (e.g., Kw an & ScovilleL 1 1976b 
ISnell et all 1 19841) . which could alter the d istribution of 12 CO in 
the region. Indeed, iMarrone et al.1 d2004l) show that broad ve- 
locity emission arises mainly from the Orion-KL region, while 
much of the narrower emission arises from the PDR excited by 
the M42 Hn region. 



5.3. Self-absorption profiles 

As noted in Sect. 13.21 all 12 CO transitions analysed in this paper 
are affected by self-absorption, which is likely to be due to cold 
gas surrounding a warmer component (Phillips et allfl98ll) . In 
particular along the infrared dark cloud (Fig. [8]), the 12 CO(6-5) 
line has blue-skewed profiles in the north-east (towards the EGO 
and the IRDC positions) and the red-skewed ones towards the 
south-west (the hot co re). Blue- and red-skewed profiles (e.g., 
Mardones et all[T997l) are commonly interpreted as due to rota- 
tion or outflow motions (which should produce equal numbers 
of red and blue profiles, and could therefore explain the profiles 
detected towards the EGO and IRDC position, the hot core and 
the Hn region, Fig. |6j or to infall (which should produce profiles 
which are skewed towards the blue, e.g. towards the EGO and 
IRDC position) or to expansion (which should produce profiles 
skewed towards the red and could be responsible for the 12 CO 
spectra of the hot core and the Hn region). From the PV diagram 
of the 12 CO(6-5) line, we do not have any evidence of global 
rotation towards the hot core and the red-skewed profile can be 
interpreted in terms of expansion or outflow motion. Similarly, 
from Fig. [10]we see that the emission does not follow a perfect 
expanding spherical shell which might imply that the rotation is 
on the origin of the self-absorbed profile seen towards the Hn 
region. Finally, the blue-skewed line profile detected toward the 
EGO and IRDC positions are typical of infall motion. 



6. Conclusions 

To study the effect of feedback from massive star forming re- 
gions in their surrounding environment, we selected the region 
G327. 3-0.6 for large scale mapping of several mid-7 12 CO and 
13 CO lines with the APEX telescope and of the high-7 13 CO(10- 
9) transition with the Herschel satellite. Our results can be sum- 
marised as follows: 
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1. Maps of all transitions are dominated by the PDR associated 
with the Hn region G327.3-0.5; mid-7 12 CO and 13 CO emis- 
sion is detected along the whole extent of the IRDC; 

2. Mid-7 transitions show rather extended emission with typi- 
cal excitation temperatures of ~ 30 K and column densities 
of some 10 21 cm" 2 ; 

3. All observed transitions are detected also in the inter-clump 
gas when averaged over large regions. The inter-clump emis- 
sion is compatible with LTE emission from a gas at 44 K, and 
with a 13 CO column density of 2xl0 15 cm" 2 , thus suggest- 
ing that the inter-clump is composed of high-density clumps 
with low filling factors; 

4. The warm gas traced by 12 and 13 CO(6-5) is only a small 
percentage (~10%) of the total gas in the infrared dark cloud, 
while it reaches values up to ~35% of the total gas in the ring 
surrounding the Hn region; 

5. The 12 CO and 13 CO ladders are qualitatively compatible 
with PDR models for high density gas (n > 10 6 cm" 3 ) and, 
in the case of 13 CO, suggest that the emission comes from a 
large number of clumps; 

6. The detection of the 13 CO(10-9) line allows to give stronger 
constraints on the physics of the gas by breaking the degen- 
eracy between density and temperature (typical of 12 CO and 
13 CO transitions) in the high temperature-low density part 
of the T - n plane. 
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